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ABSTRACT 


This is Part 11 of “A Comprehensive Report on Nineteen Condensation 


Nuclei,” which contains tabulations of the growth times of these nuclei_ 


under various assumed supersaturations. The growth time of each 
nucleus at 3 given supersaturation is ottained from a newly developed 
groeth rate equation which utilizes a single parameter to denote the 
hygroscopicity of the nucleus and which includes an additional term 
to reftect the inefficiency of the condensation process. This volume 
and its predecesscr serve as a valuabie reference fcr the experimental 
cloud physicist in his study of the growth behaviors of condensation 
nuclei and fer the field cloud nhysicist in his selection of the 
proper artificial nuclei for warm fog or cloud modification. - 


hi gg hsreceds nm cen catetee FE eee 


e 


: CONTENTS 
Pace 
Re “INTRODUCTIONS 2 2. Fy Sy PES, Be 2. BUG ten: cine, © ] 
‘GROWIH RATE EQUATIONS 2... ...0-00 2008 wean ees | 
THE LAYER CONCEPT 2... ee ee Fat ie tees 2 
_ DERIVATION OF THE EQUATION. 2. 2 2 we ee ee 3 
; Condensation: ..  . a eS eee, ym ee 3 
OGETUSION. 6.5 Se: de Be SS a We ee ee a ee ae een aw NG “4 
Thermal Conductivity ..........-266888268-6 %§ 5 
The Clapeyron-Clausius "quation.......... aids 4S 5 
7 Hygroscopicity . 2... 6. ee ee ee ee ee ee ee 6 
The New Growth Rate Equation .........+..+24.-. ? 
NUMERICAL SOLUTION 2 2 ee eae eee: 
ILLUSTRATICNS AND-COMPARISCNS . 2... ..2.-2.-0000- 9 
=" CONCLUDING REMARKS 2... 2. Present Coase Pclerenanen toa 19 
TABULATIONS OF GROWTH TIMES FOR NINETEEN CONDENSATION 
NUCLE ! : : 
|. Lithium Chloride ....... RAY ahs Nace 21 
2. Ammonium Chloride ........4.6. «Sle ee ae yee Tee 34 
3. Magnesium Chilcride .......-e- ao fect ot Sh aS Me 47 
- 4, Sodium Chloride... 6. 6 2 we ee ee ee ew ee 60 
5. Calcium Chloride ........ rise do: qa, Shee i Whe on tht 73 
6. Potassium Chloride ...... oA Say a ed GY aah ae ot te 88 
7. Ammonium Nitrate ....... oe Seat te cas bm wene tee 08 NOU 
8. Sodium Nitrate . 2... we ee eee ee we ewe 114 
9, -Sodium Bromide ......+.... divin pares, ait CRA 
10. Zinc Nitrate ... 1.1. 2s ee ew eee aes cel et ce we 2 AO 
||. Potassium Bromide... .... 2s sk, Hew es ae 1D2 
[ 12, Ammonium Sulfate ... . ee we Bs Sip ey th Sagena bs a ae SRS 
: 13. Barium Chiorid2 . ... 2. eae ot a eta to a ae ee TB 
BO 14, Sodium Sulfate .......... Dances deena Ne Hee 
15. Potassium Nitrate. ....... st. He aetca si Oe les He 4s) SLOT 
E 16, Magnesium Sulfate. ...... oh ad. Gh a tes etes, He doo ped 
_ 17. Potassium lodide . . 1 1 wee ee ee ee te ew ee) «224 
, KBs Zine Sul fates: 20 405. oo ek GL oi ee, Aa ee oes Gn eit 256 
is £9. = CUDIC-SUTOTG. Rte a ee ee gs a. Soe gtd a” SS 
| 
- LLTERATORE ‘GIleD: «ta saa SOs Ba ee Ge SS a eee ew 209 
i v 


INTRODUCT 1 ON 
This is the second part of a two-part report on the study of IY etec- 
trolytic condensation nuclei. Part | [1] deals with their equili- 
brium droplet radii versus saturation at 20C and one atmosphere. 
Utitizing the values of Fart 1, this portion extends the study to 
their growth times, an important parameter for determining their rel- 
ative efficiencies in warm fog or cloud modification. . 
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i To compute the growth rate (from which the growth time can be obtained), 
an expression *s developed on a more physically consistent basis than 
ro - heretofore presented. This expression ‘formally takes into account 

the hygroscopicity of a nucleus in terms of its water activity as wel! 
: as the condensation process in a thin hypothetical layer [2] that 
: : envelops the nucleus. The growth times of each of these 19 corden- 

“sation nuclei are presented in seven tables, corresponding to nucleus 

; mz ses between 10-!2 and 1076 grams, and each table contains six 
! : supersaturations:between 9.01% and 5%. {+ is expected that these 
tables, covering nearly all the possible combinations of nucleus sizes 
“and humidity conditions, wiil find practical application in warm fog 
or cloud seeding experimenis. 
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GROWTH RATE EQUATIONS 


: There are a number of steady-state versions of the growth rate equarion 

[3-13]. These versions are all derived from the same se: of .classicai 
transport equations, with varying degrees of mathemetical = d physical 
approximations. Most of them are used in connection with numerical 
modeling, and they serve this-purpose reasonably wel! in simulating the 
growth of a population of droplets. 


In theoretical work on the growth of a single particle, such approxi- 
mations should be avoided and the physical processes involved fully 

’ accounted for. Two processes of major importance which have not been 
handied adequately in the growth of an isolated solution droplet in 
2 Stationary medium of moist air are hygroscopicity and condensation. 
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ln the pasv, Fygroscopicity has been represented by either the van't 
Hoff factor [5] or the molal osmotic coefficient [14]. Both require 
additional parameters to form the solute term. In accord with modern 
practice in physical chemistry, the term can now be completely described 
by a single parameter; nameiy, the water activity [15]. Condensation, 
on tha other hand, has generally been lumped with diffusion of warfer 
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vapor in air to make up the so-called modified [16]* or compensated [3]** 
diffusion coefficient and has not been recognized as a distinct process 
as it should have. As a result, the physical process of condensation 
becomes intractable in the usual growth rate equation. Moreover, in 

the case of electrolytes, we suspect that how efficiently water mole- _ 
cules will condense on a liquid solution droplet may depend upon the 
_ionic property and concentration of solution at a given temperature ana 
pressure. : * 


7 THE LAYER CONCEPT 
To render the concensation process tractable mathematically and phys- 
ically, we employ the concept of boundary layer over a solution drop- 
let. Langmiur [2] used this approach in his study of hydrogen dis- 
sociation in a Light bulb. Borovikov, et al. Lil] called it the 
parietal layer and came up with a modified coefficient. Fukuta and 
Walter [17] introduced the idea of Knudsen flow in the layer and de- 
rived separate correction factors to the diffusion coefficient and 
thermal conductivity. Since the layer concept has been well accepted, 
it is proposed to treat diffusion and condensation as two distinct 
processes. : 
Outside the boundary, the transport of water molecules, being caused 
by an imbalance of partial vapor pressures between the ambient air 
and the layer, is governed by diffusion theory in gas kinetics. As 
is obvicus, this imbalance is brought about by the hygroscopic property 
of electrolytic solution. Inside the layer, to the author's know ledge, 
there is yet no adequate theory to deal with the interface problem 
ot condensation, except for the use of 3 nondimensional "fudge fac- 
tor, called condensation coefficient, to denote the inefficiency of 
the condensaticn process LI8]. Under the circumstances, the exisving 
theory of gas kinetics’will be utilized to formulate a condensation 
term to be included in the growth rate equation. 


As condensation takes place and the latent heat of condensation is 
iiberated, the droplet becomes warmer than the ambient air and so does 
tne layer due to heat conduction, both convection and radiation being 


¥ 

Rooth's formula for determining the thickness of an air layer is at 
variance with that given by Mordy [8]. The Mordy version appears to 
be correct. 


* 


¥% 
Fukuta and Walter [17] consider the use of the Cunningham factor as 
unjustified. ; 
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negligible in the present case [19]. Following the scheme. outlined 

by Langmuir [2], it is possible to derive a separate term tor neat 
transfer in the layer, but according to his calculations the tespera- 
ture at the outer boundary of the !fayer is not significantly different 
from the droplet temperature under normal atmospheric pressures. 

Thus, we may dispense with the Feat transfer term in droplet growth 
and assume that both the droplet 3nd-layer are at the same temperature. 


DERIVAT SON OF THE EQUAT!ON 


The underlying assumptions leaciig to the use of the classical taws 
of diffusion and thermal conduc~-‘vity have been discussec in some de- 
tail in the literature [19, 10, 14] and will not be repeated. We 

may now proceed to develop the new growth rate equation. 


Condensation 


By the kinetic theory of gas [.%]} on effusive tlow, we may formulate 
an expression for the amount o7 material congensed per unit time onto 
a spherical surface as follows: 


RV 
dm _ 2 ; _ 
a 4m B( 5 j ‘D4 p.)s (1) 


where m is the mass of the solution droplet, rg its radius, T its 
temnerature, + the time, Ry the gas constant of water vapor, po, the 
vapor density of the parietal 'ever whose thickness is A at Ty, pr 
the equilibrium vapor density at T, over the Groplet, and B the con- 
densation coefficient (which Ke- ner [18] called the evaporation coef- 
ficient) that denotes the fraction of impinging water molecules that 
do condense. There are some ev>erimental data on this coefficient 
for pure water [17], but none ‘or electrolytic solution. Unlike the 
diffusion coefficient, it cannct be calculated as no pertinent theory 
has been developed. For numerical computation in this paper, a value 
of 0.045 is used for no other reason than that it is given for 20C. 


The above equation of condensatio. can ve expressed in terms of the 
droplet radius 


dr RT_ 1/2 
Jey. Bua Vo = F 
dt Py. 27 (0, Pp)» (2) 


where p, is the density of the droplet. 


Sii.ots for the vapeoc density of tee layer, we cbtain 
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t* should, however, be ncted that <2} coes not hold when a saturates 
Sotution grociet has a radius comparable tc er less than the layer 
thickress; hence, 3c 


ubie Aitken mactei are excluded. Biffusiona! 
h this type cf nuclei may be handled a 
by “cisch ang Coliins L2a). 


re 
and condenseticonai ¢ 
the theory devetooed 


Siffusion 


Ustsise the laver, Fick's law of diffusien is given by 


{4} 


where r=(r, + 3) is the radius of ihe iayer, 04, the combined density 
of the layer ans drople*, D the diffusivity of water vapor in air, 
and p the vapor density at tke ampient temperature T. Stnce we deal 
only with larce and giant nuclei, it is easily seen that Tx T@ and 
XM S5s*:tution of (3) into (43 gives: 


am ca 


a (5) 


This ecuatior expresses the continuity of vapor diffus’on and conden-~ 


sation. However, the layer tocetner with the creplet, being at a 
higher temperature than the ambient due to tne release of the ‘arent 
heat of condens:ticn, resists the diffusive flow of vapor toward the 
zlayer due to the .:.calance of vapor pressures. To incorcorate this 
efiect, we divide and multiply the right-hand side of (5) by Poll ), 
tee saturation vapor density at the droplet temperature, 


y 
noe _ : dot, (ome Pa ‘ 27 ve dr P. J. 6 
—y = ay Atay BS ay - pTW7rtd-: ; 
dt p. eft Bo (Tt Rote dt ag? 


This eauation takes into account the fact that the droplet together 
with the layer experiences a slight reduction of the vapor gradient. 
{t may be observed that the last term in the bracket denotes the 

equilibrium saturation ratio over the droplet at the droplet temper- 
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ature. !f this temperature is known, (6) can be solved for the rate 
af droplet? growth. Usually, neither the d-cpiet temperature nor the 
derived quantity, 9 (T,), is know, and then ii is necessary to 
estatlish relationships between the known quantities and the unknowns. 


Therma! Conduczivity 


The law of thermal conductivity, as appiied to meist air, is given by 


dr _K . - ‘ 
"at ~ Tp, ‘'r T), 2 (7) 


where K is the thermal conductivity of water vapcr, i the latent neat 
of condensation, and T the ambient temperature. The droplet temper- 
ature can now be expressed in terms of the ambient temperature, 


rlp rLp 
7m d dr ae oO Gr 
gr oR ao st al a: (8) 


Let € = rlpg/KTidr/ct), and we have 


T. = Til . &). (9) 


The magnitude of € unoer a given supersaturation and temperature can 
be estimated. for instance, at 1% supersaturation and 20C, rdr/dt is 
approximately of the order of 0-7 in the early state of growth, L 
about 6 x 104, K about 4 x 10-7, and T about 3 x 102. Then, € is at 
feast of the orcer of 1077. However, as the droplet grows larger, & 


will become proportional ly sma! ler. : 


The Clapeyror-Clausius Equation 


To convert the saturation vapor density over a plane surface of pure 
water at the droplet temperature to that at the ambient, we make use 
of the Clapeyron-Clausius equation in the following form: 


Pea, (10) 


where V is the specific volume of water vapor, the volume occupied by 
liquid water droplets being negligible. With the aid of the perfect 
gas law, we may integrate (10) in the same manner as shown in the 

literature [7] to the saturation vapor pressure over a plane surface 


of pute water at the droplet tempe-ature, Po(T,), as 6 funct.on of 
that at the ambient, po{7). From (9) and the perfect gas law, an ex- 
press:cn in terms of seturation vapor densities is obtained. Dropping 
the oarenthesized ambient temperature notations, we have 


eee 1+ € 


0, exp a E/R TC + &)] 


(l= 5) 9 exp(L E/RT). | CH) 


Substitution of (I!) into (6) results in a term, 9/5, which is equi-- 
valent to the ambient saruration ratio, usually denoted by S.. The 
other simitar term, P-/po{T,), is the equilibrium saturation ratio 

at the droplet temperature, S,(T,), which will be dealt with in the 
next section. With the new notation (6) becomes 


Do p 2 
dr _ ° _ id 2m dr 
tas rR CS. aac €/2) (7) at 
°o Vv 
- {| - &) S(T) exp(L e/R 1) 1. : A (12) 


Hygroscopicity 


According to modern practice in physical chemistry, the proper para- 
_ meter to represent the hygroscopicity of a substance is its water 
activity, a, which varies with concentration but which remains sen- 
sibly constant under moderate temperature variations. The water 
activities of the 19 condensation nuclei are tabulated at different 
molalities in Part |! [ij of this report. The equilibrium saturation 
ratio at the droplet temperature expressed in terms of the curvature 
and solute effects is, according to Low [21], 


S = . ‘ ' 4 
S(T) arexp (20 /rp gRiT 


arexp [20'(! - &)/ro RTI 


= aeexp (20'/ro JR T)exp(-25! E/rogRT), (13) 


where of is the surface tension of the solution droplet, assumed con- 
stant with moderate temperature variations. 


The first two factors on the right-hand side of (13) are sGen te de- 
_noie the equilibrium saturation ratio at the ambient temperature, See 
' Tre above expression then becomes 


= »(-20! - 
S(T) S expt wa E/rp gRT). ; (14) 


’ The New Growth Rate Equation e 


The desired equation for the growth rate oi a single isolated solution 
droplet in a stationary medium of humid air follows upon substitution 
of (14) into (12), 


Dp o. 1/2 
dr _ "Oo a d ,27 dr . 
oF = ae {S, CI 6/2)¢ RT , dt 
d Oo ov 
L 20! ; 
-(l- £18 eo &G T - Tee (15) 


This equation now takes into account all the important physical pro- 
cesses contributing to rhe. growth of a solution droplet. Being more 
physically complete, it can be reduced to the various versions in the 
literature. Silverman and Kunkel's version [13] will result if we 
drop the (I - £/2), (I. - &), and (20'/rpgRyT) terms and use instead 
their approximate expression for the solute effect. Neiburger and 
Chien's version [9] wili be obtained it we further neglect the con- 
densation term and replace the solute effect with their consolidated 
hygroscopic term. Others will follow readily if we further expand 
the exponential factor in series and discard the higher-order terms. 
However, it wil: be well to remember the implications resulting *rom 
these approximations, since we are dealing-with rather smal] numericel 
quantities in the growth rate equation. After all, the supersatura- 
tions which are normally encountered in a natural fog or cloud are 

in the_neighborhood of 10°“. For this reason, quantities of the order 
of 10°? or greater should be retained. 


NUMERICAL SOLUTION 


An analytical solution of (15) appears te be an insurmountable problem, 
An iterative technique is used to calculate growth rates (from which 
growth times are derived), utilizing the equilibrium radii at c.ffer- 
ent saturation ratios of the various electrolytic nuclei from Part | Eta 


The coefficients and constants used in numerical integration at one 
atmosphere are listed below: : 


wD 
il 


1.72 x 107° g/cm? 


o 
R, = 4.615 x 10° ergs/g ° : 

T = 293°K 

D = 0.257 cm/sec | tO i. 
K = 4,132 x 10> cal/cem sec °K (in water vaopr) : 

L = 586.0 cal/g 

B = 0.045 


Both droplet density and surface tension are allowed to vary with con- 
centration. 


On the basis of the numericai quantities listed above, it is possible 
to reduce the algebraic complexities in numerical computation by 
making some approximations, The curvature correction term due to 
heating, 20'/rpgR,T, in the exponent is at least four orders of mag- 
nitude lower than the first term in the exponent, L/R, T, and can, there- 
fore, be dropped, The heating term, —, in (I - &), is of the order 

of 1073 during the early stage _of growth under moderate supersatura- 
tion, but will decrease to 107°, as claimed by Neiburger and Chien [9] 
and by Silverman and Kunkel [13], only when a solution droplet be- 
comes sufficiently large and, hence, rather dilute. Since we are- 
‘interested in the detailed life history of droplet growth, it would 
not be reasonable to discard a term which is only three orders of 
magnitude lower than unity. 


In contrast, the (| - €/2) factor in the condensation term may be 
approximated by unity in view of the smaller magnitudes of the other 
factors, when compared with natural supersaturations. However, it 
may be noted that the smaller the condensation coefficient is, the 
more significant the condensation term becomes. With the approxima- 
tions discussed above, (15) is rewritten as 


- (I - &) S_exp(L €/R T)] (16) 


‘eae eaNae ad 


oA ante eae ag Ee ere 
tiiinieaateitiinite "as Ninel iinet shina tia 


es es 


eas nme PE RE Te A LOE RATT ER CARI 
. . 


Sper io eee CME peer, 


which can be further manipulated to yield 


dr _ "o ied LE Dj2n , 
ao Cs, - (1 - €) S. exp Rp +5 RT C17? 


~ To facilitate numerical computation of the growth rates, it is neces- 


sary to transform (17) into a form amenable to iteration. Upon re- 
placing € with Lro4/KT(dr/dt) and letting 


a Pq Di 2a 
A= Dp, ‘" + ae 
B= hoe tt 
C= KR T7/rp jit 
u = dr/dt, 


a deceptively simple expression. fol lows 


S_- Au 
fo) 


u=C In( 


With the values of S. and r as given in Part | [Ct] and an assumed 
value of 54, Successive approximation can be applied to (18) to 

obtain u at each discrete r, from which the growth times are computed 
for the intermediate values between the siiccessive r's, the initial 

r being that of a saturated droplet for a given nucleus mass. Except 
for the initial growth time which is calculated at the initial growth 
rate for a half interval, the subsequent values are for full intervals. 
At the initia! size, the growth time is taken to be zero; this im- 
plies that it takes negligible time for a dry particle to become a 
saturated droplev. 


ILLUSTRATIONS AND COMPARISONS 


Having developed a more rigorous exoression for the growth rate of 

a liquid solution dronlet, it seems apropos to present a few examples 
to illustrate the use cf the equation and to compare the data so gen- 
erated. Two figures are presented for the illustrative purposes. 


ae 


Figure | is a study of lithium chloride (LiCl) at 1% supersaturation 
and 20C. LiCl! is illustrated in four different nucleus masses cf 
1or7'," 1079, 1o7 tl, and 107 !3 grans whose crystal sizes in microns are 
shown as black dots along the racius axis. A graph suc. as this is. 
valuable to weather modification as it enables a field cloud physicist 
to select, in a given time frame, what optimum size of nucleus to use 
under a giver humidity and updraft condition. It is quite obvious 

- that larger nuclei enjoy the advantage of being able to grow into 
larger droplets in a relatively shorter time. However, when residence 
time is not a factor, al! growth curves eventually join together one 
after the other. | . : 


- tn Figure 2, the growth times of four different nuclei of the same ; be 


mass at 1% supersaturation and 2CC are given. Their dry particles 

ure about the same size, but at 9.1 second in a supersaturated environ- 
ment their different hygroscopicities have already set them apart. 
Again, when residence time is not a factor, all nuctei, irrespective 

of their initial sizes and hygroscopicities, grow to the same size. 
With the aid of a graph such as this, it is possible to determine 

the proper condensation nucleus to use. 


The next six figures are offered for comparative purposes. Figures 

3 and 4 are comparisons with the previously published theoretica! 
values, and Figures 5 through 8 with empirical values. As-may be ex- 
pected, the new growth rate equation will give a lower rate of growth 
and, hence, a longer time, than the other versions. 


As noted earlier, (15) can be mace identical to Neiburger and Chien's an 


f9] version. Figure 3 shows such a comparison in terms cf growth 12 
rates at !% supersaturation and 5C. A NaCl particle of 2-656 x 10 
gram in mass or | micron in radius is used for this purpose. It is 
seen that the new equation predicts a lower rate, about 52%, during 
the early stage of growth and less than 1% when the droplet reaches 
100 microns. Again, the same statement also applies to a lesser 
extent to Figure 4, which is a comparison with Jiusto's [22] values 
in terms of growth times at 0.1% and |.0% supersaturations and 20C. 
The NaCl nucleus used has a mass of 107!! gram. Although the curves 
in Figure 4 do not depart from one another appreciably, about 11% on 
the average, it should be noted that such a comparison is by no 

means fair since the accumulative etfect of successive summations has 
completely masked their differences in the beginning of grow*h, 
Equation (15) approaches Fletcher's C10], which Jiusto [22] utilized 
for his computations, as the droplet grows large-. In infinite dilu- 
tion, it becomes Mason's [7] or Byers' [14]. 


In using empirical data, we have to be satisfied with a comparison 
of the growth time despite its drawba:«. Moreover, since th..e data 
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are presented in graphicat form, there is the added difficulty of 
extracting pertinent values from small graphs. Nonethele.s, compar- 
isons are made with three of Dennis' [23] six figures which illustrate 
tie growth times of submillimeter drops of sodium chloriue, ammonium 
nitrate, and calcium chloride 41 various subsaturations and with one 
of Ye's [24], which deals with sodium chloride, as the others are too 
smal! to be reproducible. 


With the exception of NaCl, ths new equation predicts a higher than 


actual rate of growth, as is rz.dily seen in these figures. To test 
whether theory generally overes:imates other electrolytic nuclei, 
Table | is presented, which sh..s a comparison tetween droplet sizes 


at a fixed length of time (Coil. in 9) and a comparison between equili- 
brium droplet sizes regardless of the length of time (Column |!) under 
the same subsaturations and tem ,eratures, as specified by Ye [24]. 
Although the percentage differzuces are somewhat disconcerting at 
first glance, it should be bors2 in mind that the inherent difficulty 
in performing this type of experiment, especially in the micron-sized 
range, is quite formidable. hcwever, some comfort may be drawn from 
the simple analysis that theor, either overestimates or underestimates 
half of the time. Near saturzeiion, experiments appear to confirm 
theory. If we can assume that their data are reasonably consistent, 

a question may be raised as to why theory, if defective, does not con- 
sistently overestimate or consistently underestimate. The answer ap- 
pears to lie in the use of the appropriate condensation coefficients 
for different electrolytic nuciei, a problem yet unexplored in chem- 
ical physics. 


CosCLUDING REMARKS 


A new expression for the growth rate of a single isolated solution 
droplet in a stationary medium of humid air has been presented, which 
takes into account the hygroscopicity of a condensation nucleus in 
terms of a single parameter in ccnsonance with the modern practice of 
physical chemistry and which includes an explicit term to reflect 

the inefficiency of the condersation process. The expression is then 
programmed to compute the growth times of 19 important condensation 
nuclei in seven different sizes under six assumed supersaturations. 


The basic properties of the 19 electrolytic nuclei, as depicted in 
Parts | and !| of this report, will fill a need for the investigators 
of applied cloud physics, not only as a handy reference to the 

nuclei considered, but also, for those who desire to use it, as an 
illustration of a meaningful method for studying otner prospective 
seecing materials. 
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